A general method for determining the masses 
of semi-invisibly decaying particles at hadron colliders 
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We present a general solution to the long standing problem of determining the masses of pair- 
produced, semi-invisibly decaying particles at hadron colliders. We define two new transverse kine- 
matic variables, Mct ± and Mct^ , which are suitable one-dimensional projections of the contrans- 
verse mass Mct- We derive analytical formulas for the boundaries of the kinematically allowed 
regions in the (Mct ± , Mct^ ) and (Mct ± , Mct) parameter planes, and introduce suitable variables 
Dot,, and Dct to measure the distance to those boundaries on an event per event basis. We show 
that the masses can be reliably extracted from the endpoint measurements of McrZ ant ^ ^cr™ ( or 
DctJ)- We illustrate our method with dilepton ti events at the LHC. 

PACS numbers: 14.80. Ly,12.60. Jv,11.80.Cr 



The ongoing run of the Large Hadron Collider (LHC) 
at CERN will finally provide the first glimpse of physics 
at the TeV scale. In large part, the excitement surround- 
ing the LHC is fueled by the anticipation of the unknown: 
no one knows for sure where or how the first signal of new 
physics beyond the standard model (BSM) will show up. 
Yet, complementary and independent arguments from 
particle physics and astrophysics suggest that the best 
place to look for new physics is a channel with missing 
transverse energy fx, caused by unseen new particles 
contributing to the dark matter of the Universe. 

Unfortunately, the study of missing energy signatures 
poses a tremendous challenge at hadron colliders like the 
LHC. The first fundamental difficulty is related to the 
very nature of hadron colliders, where in each event the 
partonic center-of-mass energy y/I and longitudinal mo- 
mentum p z of the initial state are unknown. To make 
matters worse, the lifetime of the dark matter particle 
is typically protected by a new parity symmetry, which 
guarantees that in every event the missing particles come 
in pairs, thus proliferating the number of unknown pa- 
rameters describing the final state event kinematics. 




FIG. 1: The generic event topology under consideration. All 
particles visible in the detector are clustered into three groups: 
upstream objects U with total transverse momentum Ut, and 
two composite visible particles Vi (i = 1,2), each with invari- 
ant mass mi and total transverse momentum piT. 



The generic topology of a "new physics" Ifir event is 
sketched in Fig. [T] Consider the inclusive production 
of an identical pair of new "parent" particles P. Each 
parent P decays semi-invisibly to a set Vi (i = 1,2) of 
standard model (SM) particles, which are visible in the 
detector, and a dark matter particle C (from now on 
referred to as the "child") which escapes detection. In 
general, the parent pair is accompanied by a number of 
additional "upstream" objects U (typically jets) with to- 
tal transverse momentum Ut- They may originate from 
various sources such as initial state radiation or decays of 
even heavier particles. We shall not be interested in the 
exact details of the physics responsible for U, adopting a 
fully inclusive approach to the production of the parents 
P. Given this general setup, the goal is to determine in- 
dependently the mass M p of the parent and the mass M c 
of the child in terms of U, V\ and V%. 

In the past, several approaches to this problem have 
been proposed, but each has its own limitations. For 
example, the classic method of invariant mass endpoints 
[HlH only applies when the visible SM particles in Vi arise 
from a sufficiently long decay chain. Attempts at direct 
reconstruction Q of the children momenta are again lim- 
ited to long decay chains only. In this letter, we shall con- 
sider the extreme, most challenging example where each 
visible set Vi consists of a single SM particle of fixed mass 
m, . A perfect testing ground for this scenario is provided 
by dilepton ti events (already observed at the LHC 
and we shall use that example in our numerical illustra- 
tions below. The role of the parent P (child C) will be 
played by the SM VF-boson (SM neutrino), each Vi is a 
SM lepton (e or fx), while U is composed of the two &-jets 
from the top quark decays, plus any additional QCD jets 
from initial state radiation (ISR). 

For such extremely short decay chains, the only viable 
alternative at the moment is provided by the methods 
based on the Mti variable @. There, at least in princi- 
ple, the individual masses M p and M c can be determined 
by observing a "kink" feature in the Mti endpoint as a 
function of a hypothesized trial mass M c for C [6] , or by 
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FIG. 2: Decomposition of the observed transverse momentum 
vectors from Fig. [1] in the transverse plane. 



exploring the Ut dependence of the Mt2 endpoint Q- 
Compared to those Mt2 approaches, our method here 
has two advantages. First, it is simpler - it uses only the 
observed objects U, Vi and V2 in the event and makes 
no reference to the missing particle kinematics (or mass). 
Second, it is more precise, since it utilizes the whole kine- 
matic boundary of the relevant two-dimensional distri- 
bution and not just the kinematic endpoint of its one- 
dimensional projection. We proceed in three easy steps. 

Step I. Orthogonal decomposition of the observed trans- 
verse momenta with respect to the Ut direction. The 
Tevatron and LHC collaborations currently use fixed axes 
coordinate systems to describe their data. Instead, we 
propose to rotate the coordinate system from one event 
to another, so that the transverse axes are always aligned 
with the direction Tji selected by the measured upstream 

transverse momentum vector Ut and the direction Tj_ 
orthogonal to it (see Fig. [2]). The visible transverse mo- 
mentum vectors from Fig. Q] are then decomposed as 

P* T \i = 7JT (per ■ Ut^J Ut, (1) 

PiT ± = PiT ~PiT t = Jj2^ T X (piT X U T ^j ■ (2) 

Step II. Constructing the transverse and longitudinal 
contransverse masses Mqt ± and Men ■ Our starting 
point is the original contransverse mass variable Q 



M C t = y m\ + m\ + 2 (e 1T e 2 T + Pit ■ P2t), (3) 
where eiT is the "transverse energy" of Vi 

e lT = V^ + IPd 2 - (4) 

For events with Ut = 0, Mct has an upper endpoint 
which is insensitive to the unknown Vs, providing one 
relation among M p and M c d, Q 



M^ X {U T = 0) = ^/m 2 + m 2 + 2 mi m 2 cosh + 



where 



sinh Q 



2M p m, 
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(5) 



(6) 



X(x,y,z) = x + y + z — 2xy — 2xz — 2yz . (7) 



Unfortunately, the Ut = limit is not particularly 
interesting at hadron colliders (especially for inclusive 
studies), since a significant amount of upstream Ut is 
typically generated by ISR (and other) jets. One pos- 
sible fix is to use all events, but modify the definition 
([3]) to approximately compensate for the transverse Ut 
boost Q. One then recovers a distribution whose end- 
point is still given by (0. Alternatively, one could stick 
to the original Mct variable, and simply account for the 
Ut dependence of its endpoint as 



M^ x (Ut) = yj m\ + m 2 + 2mim 2 cosh (2r, + (i + (2) 

(8) 

where Q were already defined in ([6]), and 
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(9) 



Our approach here is to utilize the one-dimensional 
projections from eqs. ()1|2|) and construct one-dimensional 
analogues of the Mct variable 



M C t ± = yml+ml + 2(e 1T± e2T ± +PiT ± -P2T ± ),{W) 
M C t 1{ = \jm\+m\ + 2 (e^ e 2 T [{ + Pit,, ■ P2T n ) , (11) 
where the corresponding "transverse energies" are 



mf + \& T± | 2 , e iT » = Jrf + \p lT ,, | 2 . (12) 



The benefit of the decomposition ()10ll 1[) is that one gets 
"two for the price of one" , i.e. two independent and com- 
plementary variables instead of the single variable ([3]). 

The variable Mct ± in particular is very useful for our 
purposes. To illustrate the basic idea, it is sufficient to 
consider the most common case, where Vi is approxi- 
mately massless (m, = 0), as the leptons in our ti ex- 
ample. A crucial property of Mct± is that its endpoint 
is independent of Ut- 



M 2 - Ml 
= — ^ — -, V U T . 



M„ 



(13) 



In fact the whole Mct± distribution is insensitive to Ut- 

dN dN 

N 0± 5(M C Tj + (N tot -N J— , (14) 
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where Nq ± is the number of events in the zero bin 
Mct ± — 0. Using phase space kinematics, we find that 
the shape of the remaining (unit-normalized) zero-bin- 
subtracted distribution is simply given by 



dN 



= —4 McT ± lnM C T ± 



dM CT± 

in terms of the unit-normalized Mct± variable 

Mct, 



M C t, 



(15) 



(16) 
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FIG. 3: Zero-bin subtracted Mct ± distribution after cuts, for 
ti dilepton events. The yellow (lower) portion is our signal, 
while the blue (upper) portion shows ti combinatorial back- 
ground with isolated leptons arising from r or b decays. 
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FIG. 4: Scatter plots of (a) Mct ± versus M C t 1{ and (b) Mct ± 
versus Mct, for a fixed representative value Ut = 75 GeV. 
The solid lines show the corresponding boundaries defined 
in ((20j) and p3|). for the correct value of Mgt± and several 



The observable Mct ± distribution for our tt example 
is shown in Fig. [3J for 10 fb" 1 of LHC data at 7 TeV. 
Events were generated with PYTHIA [l(| and processed 
with the PGS detector simulator [llj. We apply stan- 
dard background rejection cuts as follows jij: we require 
two isolated, opposite sign leptons with pit > 20 GeV, 
TO £+£- > 12 GeV, and passing a Z-veto \mg+e- — Mz\ > 
15 GeV; at least two central jets with px > 30 GeV and 
|r?| < 2.4; and a $ T cut of $ T > 30 GeV {$ T > 20 GeV) 
for events with same flavor (opposite flavor) leptons. We 
also demand at least two 6-tagged jets, assuming a flat 
6-tagging efficiency of 60%. With those cuts, the SM 
background from other processes is negligible [4|. 

Fig. [3] demonstrates that the Mct ± endpoint can be 
measured quite well. Since the theoretically predicted 
shape (fl"5|) is distorted by the cuts, we use a linear slope 
with Gaussian smearing, and fit for the endpoint and the 
resolution parameter. We find M™^ — 80.9 GeV (com- 
pare to the true value M™^ x = 80.4 GeV), which gives 
one constraint (|13l) among M p and M c . At this point, 
a second, independent constraint can in principle be ob- 
tained from an analogous measurement of the Mq^ x end- 
point ([8]) at a fixed value of Ut (resulting in loss in statis- 
tics), after which the two masses can be found from 



M p 



U T M^ X (U T )M, 
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(17) 
(18) 



However, the orthogonal decomposition (|10lllj) offers an- 
other approach, which we pursue in the last step. 

Step III. Fitting to kinematic boundary lines. It is 
known that two-dimensional correlation plots reveal a 
lot more information than one-dimensional projected his- 
tograms @, [HI . To this end, consider the scatter plot of 
Mqt ± vs Mct» in Fig- IUa), where for illustration we 



used 10,000 events at the parton level. For a given value 
of Mct± j the allowed values of Mct, are bounded by 



Af <&> (M CT± ) < Mct, < (M CT± ), (19) 
where (M CT± ) = and 

Mct] (Mct ± ) = M%%_ (jl — M CT± cosh r, + sinh r, 



(20) 

Fig. |4ja) reveals that the endpoint Mq^ x of the one- 
dimensional Mct, distribution is obtained at Mct ± = 

= M<$> (0) = M^l (cosh n + sinh r,) 

= {^ + U T + U T ). (21) 

Notice that events in the zero bins Mqt ± — and 
Mct, = fall on one of the axes and cannot be dis- 
tinguished on the plot. 

Now consider the scatter plot of Mct± vs Mct shown 
in Fig. IDJb). Mct is similarly bounded by 

M<£? {M CT± ) < M CT < M^ (M CT± ) , (22) 

where this time M c l j! (Mct ± ) — Mct ± and 



(M CT± ) = M^ x cosh r] + Jl- M CT sinh j 



(23) 

We see that the endpoint M™^ x of the one-dimensional 
Mct distribution is also obtained for Mct± = 0: 



Mct x = M { c h ^(0) = Af^(cosh»7 + sinh?]) = Mgff 



(24) 
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FIG. 5: Dct distributions for four different values of M p (and 
M c given from (|18p ). The yellow (light shaded) histograms 
use only events in the zero bin Mct ± = 0. The red solid lines 
show linear binned maximum likelihood fits. 
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FIG. 6: Fitted values of D^ n as a function of M p . 



Fig. |4] reveals a conceptual problem with one- 



dimensional projections. While all points in the vicin- 
ity of the boundary lines (|20|) and (|23|) are sensitive 
to the masses, the M^ x endpoint is extracted mostly 
from events with M CT± ~ Mg^, while the and 
M(jf x endpoints are extracted mostly from the events 
with Mct ± ~ 0. The events near the boundary, but with 
intermediate values of Mct ± , will not enter efficiently ei- 
ther one of these endpoint determinations. 

So how can one do better, given the knowledge of the 
boundary line ([23)) ? In the spirit of [13| . we define the 
signed distance to the corresponding boundary, e.g. 

D CT (M P , M c ) = M%j> (M CT± , U T , M p , M e ) - M CT 

and similarly for Dct« ■ The key property of this variable 
is that for the correct values of M p and M c , its lower 
endpoint D™™ is exactly zero (see Fig. Efb)): 

D^ n (Mp,M c )=0. (25) 
In that case the boundary line provides a perfectly snug 
fit to the scatter plot — notice the green boundary line 
marked "80" in Fig.@|b). While in general eq. (|2"5"j) rep- 
resents a two-dimensional fit to M p and M c , in practice 
one can already use the measurement to reduce 

the problem to a single degree of freedom, e.g. the par- 
ent mass M p , as presented in Figs. 2] and [SJ We see that 
the correct parent mass M p = 80 GeV provides a perfect 
envelope, for which D™™ = 0. If, on the other hand, M p 
is too low, a gap develops between the outlying points 
in the scatter plot and their expected boundary, which 
results in D™™ > 0. Conversely, if M p is too high, some 
of the outlying points from the scatter plot fall outside 
the boundary and have Dct < 0, leading to D™™ < 0, 
as seen in Fig. EIc,d). The resulting fit for Z?™™ as a 
function of M p from our PGS data sample is shown in 
Fig. [6j which suggests that a W mass measurement at 
the level of a few percent might be viable. 
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